home *** CD-ROM | disk | FTP | other *** search
/ Shareware Grab Bag / Shareware Grab Bag.iso / 002 / pc_pad.arc / TRANS.BAS (.txt) < prev    next >
Encoding:
GW-BASIC  |  1983-07-03  |  5.1 KB  |  102 lines
  1. 10  REM - This is TRANS
  2. 20  KEY OFF: CLS: COLOR 0,7,0
  3. 30  PRINT " Transcendental Functions of Complex Numbers "
  4. 40  REM - By Robert Braunwart
  5. 50  COLOR 7,0,0
  6. 60  DEF FNS(X)=(EXP(X)-EXP(-X))/2
  7. 70  DEF FNC(X)=(EXP(X)+EXP(-X))/2
  8. 80  PI=3.14159
  9. 90  PRINT
  10. 100  PRINT "Input Z (X,Y)";
  11. 110  INPUT X,Y
  12. 120  PRINT
  13. 130  SR=SIN(X)*FNC(Y): SI=COS(X)*FNS(Y)
  14. 140  CR=COS(X)*FNC(Y): CI=(-1)*(SIN(X))*(FNS(Y))
  15. 150  SHR=FNS(X)*COS(Y): SHI=FNC(X)*SIN(Y)
  16. 160  CHR=FNC(X)*COS(Y): CHI=FNS(X)*SIN(Y)
  17. 170  PRINT "   exp  Z =";
  18. 180  PRINT USING "+######.#####";EXP(X)*COS(Y),EXP(X)*SIN(Y);
  19. 190  PRINT " i"
  20. 200  PRINT "   sin  Z =";
  21. 210  PRINT USING "+######.#####";SR,SI;
  22. 220  PRINT " i"
  23. 230  PRINT "   cos  Z =";
  24. 240  PRINT USING "+######.#####";CR,CI;
  25. 250  PRINT " i"
  26. 260  PRINT "   tan  Z =";
  27. 270  PRINT USING "+######.#####";(SR*CR+SI*CI)/(CR*CR+CI*CI),(CR*SI-SR*CI)/(CR*CR+CI*CI);
  28. 280  PRINT " i"
  29. 290  PRINT "   sinh Z =";
  30. 300  PRINT USING "+######.#####";SHR,SHI;
  31. 310  PRINT " i"
  32. 320  PRINT "   cosh Z =";
  33. 330  PRINT USING "+######.#####";CHR,CHI;
  34. 340  PRINT " i"
  35. 350  PRINT "   tanh Z =";
  36. 360  PRINT USING "+######.#####";(SHR*CHR+SHI*CHI)/(CHR*CHR+CHI*CHI),(CHR*SHI-SHR*CHI)/(CHR*CHR+CHI*CHI);
  37. 370  PRINT " i"
  38. 380  PRINT
  39. 390  IF SGN(Y)=0 AND SGN(X)=0 THEN PRINT "   log  Z =     -/": GOTO 460
  40. 400  PRINT "   log  Z =";
  41. 410  IF SGN(X)=0 THEN T=SGN(Y)*PI/2: GOTO 440
  42. 420  T=ATN(Y/X)
  43. 430  IF SGN(X)=-1 THEN T=T+PI
  44. 440  PRINT USING "+######.#####";LOG(SQR(X*X+Y*Y)),T;
  45. 450  PRINT " i + 2k<0xE3!>i"
  46. 460  IF SQR(X*X+Y*Y)>0.75 GOTO 950
  47. 470  X1=(X)+(X^3-3*X*Y^2)/6+3*(X^5-10*X^3*Y^2+5*X*Y^4)/40+5*(X^7-21*X^5*Y^2+35*X^3*Y^4-7*X*Y^6)/112+35*(X^9-36*X^7*Y^2+126*X^5*Y^4-84*X^3*Y^6+9*X*Y^8)/1152
  48. 480  X1=X1+63*(X^11-55*X^9*Y^2+330*X^7*Y^4-462*X^5*Y^6+165*X^3*Y^8-11*X*Y^10)/2816
  49. 490  X1=X1+231*(X^13-78*X^11*Y^2+715*X^9*Y^4-1716*X^7*Y^6+1287*X^5*Y^8-286*X^3*Y^10+13*X*Y^12)/13312
  50. 500  X1=X1+1001*(X^15-105*X^13*Y^2+1365*X^11*Y^4-5005*X^9*Y^6+6435*X^7*Y^8-3003*X^5*Y^10+455*X^3*Y^12-15*X*Y^14)/71680
  51. 510  Y1=(Y)+(3*X^2*Y-Y^3)/6+3*(5*X^4*Y-10*X^2*Y^3+Y^5)/40+5*(7*X^6*Y-35*X^4*Y^3+21*X^2*Y^5-Y^7)/112+35*(9*X^8*Y-84*X^6*Y^3+126*X^4*Y^5-36*X^2*Y^7+Y^9)/1152
  52. 520  Y1=Y1+63*(11*X^10*Y-165*X^8*Y^3+462*X^6*Y^5-330*X^4*Y^7+55*X^2*Y^9-Y^11)/2816
  53. 530  Y1=Y1+231*(13*X^12*Y-286*X^10*Y^3+1287*X^8*Y^5-1716*X^6*Y^7+715*X^4*Y^9-78*X^2*Y^11+Y^13)/13312
  54. 540  Y1=Y1+1001*(15*X^14*Y-455*X^12*Y^3+3003*X^10*Y^5-6435*X^8*Y^7+5005*X^6*Y^9-1365*X^4*Y^11+105*X^2*Y^13-Y^15)/71680
  55. 550  PRINT "arcsin  Z =";
  56. 560  PRINT USING "+######.#####";X1,Y1;
  57. 570  PRINT " i + 2k<0xE3!>"
  58. 580  PRINT "arccos  Z =";
  59. 590  PRINT USING "+######.#####";(PI/2)-X1,-1*Y1;
  60. 600  PRINT " i + 2k<0xE3!>"
  61. 610  X2=(X)-(X^3-3*X*Y^2)/3+(X^5-10*X^3*Y^2+5*X*Y^4)/5-(X^7-21*X^5*Y^2+35*X^3*Y^4-7*X*Y^6)/7+(X^9-36*X^7*Y^2+126*X^5*Y^4-84*X^3*Y^6+9*X*Y^8)/9
  62. 620  X2=X2-(X^11-55*X^9*Y^2+330*X^7*Y^4-462*X^5*Y^6+165*X^3*Y^8-11*X*Y^10)/11+(X^13-78*X^11*Y^2+715*X^9*Y^4-1716*X^7*Y^6+1287*X^5*Y^8-286*X^3*Y^10+13*X*Y^12)/13
  63. 630  X2=X2-(X^15-105*X^13*Y^2+1365*X^11*Y^4-5005*X^9*Y^6+6435*X^7*Y^8-3003*X^5*Y^10+455*X^3*Y^12-15*X*Y^14)/15
  64. 640  Y2=(Y)-(3*X^2*Y-Y^3)/3+(5*X^4*Y-10*X^2*Y^3+Y^5)/5-(7*X^6*Y-35*X^4*Y^3+21*X^2*Y^5-Y^7)/7+(9*X^8*Y-84*X^6*Y^3+126*X^4*Y^5-36*X^2*Y^7+Y^9)/9
  65. 650  Y2=Y2-(11*X^10*Y-165*X^8*Y^3+462*X^6*Y^5-330*X^4*Y^7+55*X^2*Y^9-Y^11)/11+(13*X^12*Y-286*X^10*Y^3+1287*X^8*Y^5-1726*X^6*Y^7+715*X^4*Y^9-78*X^2*Y^11+Y^13)/13
  66. 660  Y2=Y2-(15*X^14*Y-455*X^12*Y^3+3003*X^10*Y^5-6435*X^8*Y^7+5005*X^6*Y^9-1365*X^4*Y^11+105*X^2*Y^13-Y^15)/15
  67. 670  PRINT "arctan  Z =";
  68. 680  PRINT USING "+######.#####";X2,Y2;
  69. 690  PRINT " i + 2k<0xE3!>"
  70. 700  X3=(X)-(X^3-3*X*Y^2)/6+3*(X^5-10*X^3*Y^2+5*X*Y^4)/40-5*(X^7-21*X^5*Y^2+35*X^3*Y^4-7*X*Y^6)/112+35*(X^9-36*X^7*Y^2+126*X^5*Y^4-84*X^3*Y^6+9*X*Y^8)/1152
  71. 710  X3=X3-63*(X^11-55*X^9*Y^2+330*X^7*Y^4-462*X^5*Y^6+165*X^3*Y^8-11*X*Y^10)/2816
  72. 720  X3=X3+231*(X^13-78*X^11*Y^2+715*X^9*Y^4-1716*X^7*Y^6+1287*X^5*Y^8-286*X^3*Y^10+13*X*Y^12)/13312
  73. 730  X3=X3-1001*(X^15-105*X^13*Y^2+1365*X^11*Y^4-5005*X^9*Y^6+6435*X^7*Y^8-3003*X^5*Y^10+455*X^3*Y^12-15*X*Y^14)/71680
  74. 740  Y3=(Y)-(3*X^2*Y-Y^3)/6+3*(5*X^4*Y-10*X^2*Y^3+Y^5)/40-5*(7*X^6*Y-35*X^4*Y^3+21*X^2*Y^5-Y^7)/112+35*(9*X^8*Y-84*X^6*Y^3+126*X^4*Y^5-36*X^2*Y^7+Y^9)/1152
  75. 750  Y3=Y3-63*(11*X^10*Y-165*X^8*Y^3+462*X^6*Y^5-330*X^4*Y^7+55*X^2*Y^9-Y^11)/2816
  76. 760  Y3=Y3+231*(13*X^12*Y-286*X^10*Y^3+1287*X^8*Y^5-1726*X^6*Y^7+715*X^4*Y^9-78*X^2*Y^11+Y^13)/13312
  77. 770  Y3=Y3-1001*(15*X^14*Y-455*X^12*Y^3+3003*X^10*Y^5-6435*X^8*Y^7+5005*X^6*Y^9-1365*X^4*Y^11+105*X^2*Y^13-Y^15)/71680
  78. 780  PRINT "arcsinh Z =";
  79. 790  PRINT USING "+######.#####";X3,Y3;
  80. 800  PRINT " i + 2k<0xE3!>i"
  81. 810  PRINT "arccosh Z =";
  82. 820  PRINT USING "+######.#####";Y1,(PI/2)-X1;
  83. 830  PRINT " i + 2k<0xE3!>i"
  84. 840  X4=(X)+(X^3-3*X*Y^2)/3+(X^5-10*X^3*Y^2+5*X*Y^4)/5+(X^7-21*X^5*Y^2+35*X^3*Y^4-7*X*Y^6)/7+(X^9-36*X^7*Y^2+126*X^5*Y^4-84*X^3*Y^6+9*X*Y^8)/9
  85. 850  X4=X4+(X^11-55*X^9*Y^2+330*X^7*Y^4-462*X^5*Y^6+165*X^3*Y^8-11*X*Y^10)/11
  86. 860  X4=X4+(X^13-78*X^11*Y^2+715*X^9*Y^4-1716*X^7*Y^6+1287*X^5*Y^8-286*X^3*Y^10+13*X*Y^12)/13
  87. 870  X2=X2+(X^15-105*X^13*Y^2+1365*X^11*Y^4-5005*X^9*Y^6+6435*X^7*Y^8-3003*X^5*Y^10+455*X^3*Y^12-15*X*Y^14)/15
  88. 880  Y4=(Y)+(3*X^2*Y-Y^3)/3+(5*X^4*Y-10*X^2*Y^3+Y^5)/5+(7*X^6*Y-35*X^4*Y^3+21*X^2*Y^5-Y^7)/7+(9*X^8*Y-84*X^6*Y^3+126*X^4*Y^5-36*X^2*Y^7+Y^9)/9
  89. 890  Y4=Y4+(11*X^10*Y-165*X^8*Y^3+462*X^6*Y^5-330*X^4*Y^7+55*X^2*Y^9-Y^11)/11
  90. 900  Y4=Y4+(13*X^12*Y-286*X^10*Y^3+1287*X^8*Y^5-1726*X^6*Y^7+715*X^4*Y^9-78*X^2*Y^11+Y^13)/13
  91. 910  Y2=Y2+(15*X^14*Y-455*X^12*Y^3+3003*X^10*Y^5-6435*X^8*Y^7+5005*X^6*Y^9-1365*X^4*Y^11+105*X^2*Y^13-Y^15)/15
  92. 920  PRINT "arctanh Z =";
  93. 930  PRINT USING "+######.#####";X4,Y4;
  94. 940  PRINT " i + 2k<0xE3!>i"
  95. 950  PRINT
  96. 960  PRINT "More data";
  97. 970  INPUT A$
  98. 980  IF A$="NO" OR A$="No" OR A$= "no" OR A$="N" OR A$="n" GOTO 1000
  99. 990  GOTO 90
  100. 1000  RUN "MATH
  101. 1010  END
  102.